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Abstract 

We discuss the counting statistics and current noise associated with the 
double Josephson quasiparticle resonance point in a superconducting single 
electron transistor. The counting statistics are in general phase-dependent, 
despite the fact that the average current has no dependence on phase. 
Focusing on parameter regimes where the counting statistics have no phase- 
dependence, we use a general relation first derived by MacDonald in 1948 
to obtain the full frequency-dependent shot noise directly from the counting 
statistics, without any further approximations. We comment on problems 
posed by the phase-dependence of the counting statistics for the finite- 
frequency noise. 

1. Introduction 

Resonant Cooper-pair tunneling, also known as Josephson quasiparticle 
tunneling, refers to transport cycles in superconducting single-electron tran- 
sistors (SSET's) which involve the transfer of both Cooper pairs and quasi- 
particles ^[21 EH- They have recently been the subject of renewed attention, 
both because of their unusual noise properties [UE3E1 an d because of their 
utility in measuring the state of a charge superconducting qubit [7J|S] . In 
terms of noise properties, it has been shown that charge fluctuations asso- 
ciated with these processes can induce a population inversion in a coupled 
two-level system (i.e in terms of its charge noise, the transistor effectively 
has a negative temperature) jHUHj- The shot-noise in the current through the 
transistor was also found to have remarkable properties By tuning the 
strength of the Cooper pair tunneling relative to the quasiparticle tunneling, 
one could effectively tune the Fano factor determining the zero-frequency 
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shot noise. It was possible to achieve a Fano factor greater than one, which 
was interpreted as a consequence of the effective charge associated with the 
transport cycle being greater than one. Perhaps more surprisingly, it was 
possible to reduce the Fano factor below 1/2, behaviour that was not fully 
explained. The finite frequency current noise also showed interesting be- 
haviour [3]- in the regime where the Cooper-pair tunneling dominated the 
quasiparticle tunneling, a coherent peak in the current noise was predicted 
at the Josephson energy. 

In this paper, we look specifically at the statistics of transferred charge 
and current noise associated with the double Josephson quasiparticle (DJQP) 
process in a SSET, a process which allows for near quantum- limited mea- 
surement (Hj, and which has been used in a recent experiment We pay 
careful consideration to the fact that the breaking of gauge invariance by 
superconductivity can make the interpretation of counting statistics more 
subtle than in the normal-metal case |1U1 1121 I13| . One obtains phase- 
dependent counting statistics here despite the fact that transport is in- 
sensitive to the overall phase difference between the two superconducting 
reservoirs. We identify two regimes in which there is no ambiguity in defin- 
ing counting statistics, and proceed to use them to directly calculate the 
frequency-dependent shot noise, without any further approximations. This 
is in contrast to previous treatments both of normal and superconducting 
single electron transistors |141 1151 11] , where an additional Markov approx- 
imation (beyond what is needed for the counting statistics) is made to 
obtain the frequency-dependent current noise. The direct connection be- 
tween counting statistics and frequency-dependent current noise we employ 
here allows one to give a simple interpretation of the latter. In particular, 
coherent peaks in the shot noise appearing at the Josephson energy can 
be directly tied to the time dependence of the counting statistics. We also 
suggest that the suppression of the Fano factor below 1/2 is related to the 
coherence of Cooper-pair tunneling. 

2. Description and Model for Josephson Quasiparticle Tunneling 

A SSET consists of a superconducting, Coulomb-blockaded island which is 
attached via tunnel junctions to two superconducting electrodes (Fig. P). 
The SSET Hamiltonian H s = H K + H c + H v + H T has a term H K de- 
scribing the kinetic energy of source, drain and central island electrons (i.e. 
three bulk BCS Hamiltonians), a term Hy which describes the work done 
by the voltage sources, and a tunneling term Ht- The charging interaction 
is He = Ec(n — AT) 2 ) where Eq is the SSET capacitive charging energy, n 
is the number of electrons on the central island, and N is the dimensionless 
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Figure 1. Schematic showing the four steps of the double Josephson quasiparticle process 
which can occur in a superconducting single-electron transistor. Circles represent the 
central island of the SSET, while the rectangles are the electrodes. Numbers indicate the 
charge of the SSET island. 



gate voltage applied to the island. We consider a SSET with identical tunnel 
junctions, whose dimensionless conductance g satisfies g/(2n) <C 1. 

The DJQP process occurs when both the transistor gate voltage M and 
drain-source voltage 2Vds are tuned such that two Cooper-pair tunneling 
transitions, one in each junction, are resonant. We label these transitions 
as n = —* 2 in the left junction, and n = 1 — ► —1 in the right junction 
(see Fig. ^). The double resonance point for these transitions occurs at 
eVos = Ecs an d Ms = 1/2. In addition, Ec/A- (where A is the super- 
conducting gap of the SSET) must be chosen so that the quasiparticle 
transitions linking the two Cooper pair resonances are energetically allowed 
(i.e. n = 2 — > 1 and n = — 1 — ► 0), whereas unwanted transitions which 
would interrupt the cycle (i.e. n = — ► 1) are not. We take Ec = A to 
satisfy these conditions; this corresponds to the experiment of Ref. |Hj. 

Assuming that all the above conditions are met, and that T <C Ec, 
transport through the SSET will be dominated by the DJQP process sketched 
in Fig.^ The standard theoretical description of this process is obtained 
by looking at the dynamics of the reduced density matrix p describing 
the charge n of the central SSET island. Here, only four charge states are 
important; further, off-diagonal terms need only be retained between states 
involved in Josephson tunneling. Letting pu t j\ = (n = i\p\n = j), we can 
represent the non-zero elements of this reduced density matrix as an 8 
component vector: 

P = ( P(-l -1) P(0,o) P(l,l) P(2,2) P(-l,l) P(l ,-l) P(0,2) P(2,o) ) (1) 

By starting with the von Neumann equation for the evolution of the full 
density matrix, treating the tunneling Hamiltonian to lowest order in per- 
turbation theory, and then tracing out the fermion degrees of freedom (see, 
e.g. Ref. |15j). one obtains the an evolution equation for p. Taking both 
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junctions to be on resonance, we have: 



j t p{t) = U^MUp(t) 



(2) 



M 



-r 











■Ej 

- % ~2 


4 2 








\ 


r 

















' 2 


1 2 













r 


4 2 


4 2 



















-r 








4 2 


4 2 














-r/2 













1 2 





■% 








-r/2 


















-ih 

4 2 








-r/2 










F 2 
L 2 














-r/2 


/ 



(3) 



U = diag ( 1, 1, 1, 1, e 



(4) 



= gA/8 is the standard Ambegaokar-Baratoff value for the Josephson 
energy emerging from perturbation theory in the tunneling Hamiltonian; 
4>L and (fin are the phases of the superconducting reservoirs. T is the rate 
associated with the two quasiparticle transitions occurring the DJQP cycle, 
and is given by the usual expression for quasiparticle tunneling between two 
superconductors f° r simplicity, we assume the two transitions to have 
equal rates. We have also made a Markov approximation to obtain Eq. ©, 
which is valid as long we do not probe the evolution of p on timescales 
smaller than h/Ec- 

The zero eigenvector po of the evolution matrix M in Eq. © corresponds 
to the stationary value of the reduced density matrix describing the island 
charge; note that in general, it will have non-zero off-diagonal elements. 
The diagonal elements of po (which are probabilities) can be obtained from 
a classical rate equation; one eliminates the off-diagonal elements from Eq. 
P|l by expressing them in terms of the diagonal elements. Letting p be the 
4-vector of these probabilities, we are lead to the rate equation: 
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p(t) = o 



(5) 



where 7 = Ej/T represents a "rate" for Cooper-pair tunneling. Eq. 
represents an incoherent model for the DJQP process, and is sufficient 
for many purposes (e.g., calculating the average current (I) = 37r/(47 + 
2r)). The full stationary value of the density matrix (including off-diagonal 
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terms) is given by: 
1 



Po 



4 7 + 2r 



7,7 + r,7 + r, 7 , iEje-Wx , -iEje^ R , -iEje^ , iEje 



-i<t>L 



o l <PL 



(6) 



3. Counting Statistics 



In order to obtain information about the number of electrons that have 
tunneled through one of the SSET junctions, the density matrix approach 
of the previous section must be embellished. The standard approach for 
single electron transistors |15l H] is to employ a "counter" scheme where one 
explicitly tracks the dynamics of, say, itil, the number of electrons which 
have tunneled through the left junction of the SSET. The Hilbert space of 
the SSET is expanded so that each state is now also labelled by a value of 
rriL, and the tunnel Hamiltonian is modified so that it raises or lowers this 
index as appropriate (e.g., by introducing an auxiliary raising operator Ft 
such that [mL,F^] = F>). The resulting equation for the reduced density 
matrix describing both n and mi is simplified if one Fourier transforms in 
the latter variable. One finds that ^.p{k) = M (k)U p(k) , where we have 
defined: 
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if n\ = 0, ri2 = 2, 
if m = 2,n 2 = 0, (7) 
otherwise. 
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(8) 

The probability distribution for rriL is thus easily obtained from the char- 
acteristic function p(k, t): 



pirn L . /) / ^ ( 



dk 

- ik, "",U:l)=. /^e^^ M (M) (9) 



p(k,t) = e M ^ t p\k,t = 0) 



(10) 



DJQPNoise2.tex; 2/02/2008; 0:10; p. 5 



6 



It only remains to specify the initial state of the transistor (i.e. p(k, t = 0)). 
The natural choice is to place the transistor n degree of freedom in its 
stationary state, and choose uil = with probability one (i.e. p(k,t = 
0) = po) ^E]- p(rriL,t) would then represent the stationary probability that 
mi electrons are transferred through the left junction in a time t. While 
this procedure is perfectly well defined for non-superconducting systems 
[TB] . it is problematic in the superconducting case. The stationary state of 
the SSET has coherence between different charge states- thus, one cannot 
have both that mi = and the coherence required in the stationary state. 
A possible remedy for this problem would be to suppress the off-diagonal 
elements of the initial density matrix; of course, this procedure needs to 
be justified. Note that given the linearity of Eq. (jlOj) in the initial density 
matrix, we can separate incoherent and coherent contributions to p(rriL,t) 
(i.e. those arising, respectively, from the diagonal and off-diagonal elements 
of the initial density matrix). 

To gain further insight, it is instructive to consider an alternate to 
the "counter" scheme discussed above. One can also obtain a tunnelled 
charge distribution for our system by following the general prescription 
for obtaining full counting statistics discussed in Ref. |3 EH; we denote 
this distribution as p'(rriL,t). The procedure here involves considering the 
effects of coupling an ideal measurer of current to the transistor. In practice, 
a gauge transformation is made to eliminate the interaction with the mea- 
surement device, resulting in a phase e lX ^ being attached to the tunnel 
matrix element t Using the fact that the perturbative calculation 

of p in the tunneling Hamiltonian has a Keldysh structure (see e.g. Ref. 
15 J, one finds that X(t) = 5 + k on the forward Keldysh contour, and 
5 — k on the backwards Keldysh contour. 5 is a phase which depends on 
the initial state of the detector; for non-superconducting systems, gauge 
invariance ensures that it does not play a role. This procedure results in 
an expression for p'(k, t) which is identical in form to Eq. (jlOj) . except that 
now the initial state is uniquely specified as po, and the evolution matrix 
M undergoes a 5-dependent unitary transformation which only affects the 
coherent contribution to p(k,t): 



Thus, one has in general phase-dependent counting statistics for our system, 
meaning that there is a dependence on the phase shift 5 introduced by 
the measurement. This phase-dependence exists despite the fact that the 
average current is not sensitive to the phase difference between the two 
superconducting reservoirs. As discussed in Ref. E3E3, the results of any 
charge-counting experiment will now be dependent on the initial state of 
the detector. We thus see that the ambiguity in defining the initial state of 
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the transistor in the "counter" scheme translates here to needing to know 
the initial state of the detector. Note that the incoherent part of p(k, t) is 
the same in the two schemes, while the coherent part differs. The process of 
ignoring the initial coherence of the density matrix in the "counter" scheme 
corresponds here to assuming an initial detector state which is completely 
delocalized in the phase 8. 

Despite the above caveats, we can determine the counting statistics un- 
ambiguously in two limits where the off-diagonal elements of the stationary 
state po are suppressed- either F S> Ej or Ej 3> T (c.f. Eq. ©). In 
these limits the coherent contribution to p(k, t) vanishes, implying that the 
counting statistics are identical in both schemes, and have no dependence 
on the phase shift 5. Note that in general, the magnitude of the off-diagonal 
elements of p in the stationary state are driven both by the size of Ej/T 
and by the population asymmetry of the two charge states involved (i.e. 
P02 oc (poo — P22))- In the limit of large Ej, it is the lack of population 
asymmetry which suppresses the off-diagonal elements. In this limit, the 
symmetric and anti-symmetric superpositions of charge states are equally 
populated, leading to a vanishing of the off-diagonal matrix elements. 



3.1. THE LIMIT Ej < T 



In this limit, the SSET effectively gets stuck in the states n = and n = 1 
waiting for the relatively slow Cooper-pair transitions to occur; we can 
think of the rate 7 (c.f. Eq. ©) as describing effective transitions between 
n = and n = 1 (i.e. 7 describes both a Cooper-pair event and the 
subsequent quasiparticle event). We find for the probability distribution 
p(m,t) =p{mi,t) = p(mR,t) (i.e. the junctions are identical) : 



p(m, t) 



ifm = 3l 



e 



-2^tSW if m = 3Z + l or m = 3/ + 2 ( 12 ) 
otherwise. 



We have neglected terms which are small as (Ej/T) 2 , and chosen m = 
at t = 0. Again, as the coherent contribution to p(m, t) vanishes to leading 
order, there is no ambiguity in defining the counting statistics. The statistics 
can be given a simple interpretation. Each Cooper-pair plus quasiparticle 
event is described by the rate 7 and a Poisson distribution. After an even 
number of 7 transitions, a multiple of 3 electrons must have been transferred 
through a junction. However, after an odd number of 7 transitions, there 
is an equal probability of having had an extra two electrons transferred 
(i.e. if we are looking at the left junction, n = — » 2 — > 1) or of having 
a single extra electron transferred (i.e. n = 1 — > — 1 — ► for the left 
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junction). Note that there is an asymmetry here which favours values of 
m which are an integer multiple of 3; this is the analog of the even-odd 
asymmetry found for the single JQP process Also, note that in this 
limit the counting statistics are identical if we use the incoherent model of 
Cooper-pair tunneling described by Eq. ®; coherence plays no role. 



3.2. THE LIMIT Ej > T 

In this limit, we can think of there being long periods of coherent Josephson 
oscillations which are interrupted by infrequent quasiparticle transitions; 
note that we still assume Ej <C Ec = A, so that it does not effect the 
energetics of quasiparticle tunneling. Defining the Poisson distribution A(Z) 
as: 

A(Z) = e- R / 2 ^ (13) 
we find up to terms of order Ej/T: 



< i 



p(m, t) 



3 (A(2Z) (3 + ^g^j + A(2Z + 1) + A(2l - 1) J if m = 31, 

\ (A(2Z + 1) + ±A(2Z + 2) (l - c ^f)) if m = 31 + l [u) 

\ (A(2Z + 1) + |A(2Z) (l - ifm = 3Z + 2, 

otherwise 



where again, we have assumed an initial state with m = without any 
ambiguity, as the coherent contribution to p(m, t) is negligible. We can 
again think of the quasiparticle transitions as effectively being described 
by a rate 7 = T/2 and a Poisson distribution. Now, however, after an 
even number of transitions (say 2Z), it is not certain that 3Z electrons will 
have been transferred through a junction; rather, due to weakly damped 
Josephson oscillations, there is also a probability to find 3Z ± 2 electrons 
transferred. Even more remarkably, there are signatures of coherent oscil- 
lations in p(m,t), this despite the fact that the initial coherence of the 
density matrix is irrelevant. As we will see in the next section, the time- 
dependence of the counting statistics has a direct impact on the frequency 
dependence of the current noise. Note that the oscillation serves to modu- 
late the asymmetry which favours m being a multiple of three. For Ft 3> 1, 
the probability of m being a multiple of 3 is given by: 

P (m = 31) -> 5 + COS o {Ejt) (15) 
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4. Current Noise and Counting Statistics 



In this section, we calculate the frequency dependent current noise Sj(to) 
associated with the DJQP process. As is standard, the capacitance of the 
junctions Cl and Cr may be used to connect this quantity (which in- 
cludes the effects of displacement currents) to the noise associated with the 
tunneling currents in each contact |14j : 

S>M = ^S ll (lo) + ^S rr (lo) + %^ (S M H + Srl(u})) (16) 



Cs Cs Cs 



where 



•„{, = lim lim f T dte iut (({I a (t + t ),Ip(t )})). (17) 



The order of limits here ensures that we are taking averages with respect to 
the stationary density matrix of the SSET; note also that we are calculating 
the classical part of the current noise, which is a symmetric function of 
frequency. Previous calculations of this quantity for both normal and super- 
conducting transistors have effectively made use of the so-called quantum 
regression theorem, which involves making an additional Markov approx- 
imation beyond that necessary to calculate the reduced density matrix p 
[bH IT5| |3] . It turns out that this is not necessary; one can directly connect 
the particle-current noise correlators to the counting statistics using a for- 
mula first derived by MacDonald |T7] . Using the definition I a (t) = dtm a (t), 
one can show: 

S a p(u) = 2e 2 oj J dt {smut) — ((m a (t)mp(t))) , (18) 

where at t = 0, mi, = itir = with certainty, and the system is described by 
the stationary density matrix pq. The integral in this expression should be 
interpreted within the theory of distributions, i.e. J* dt sinuit = We 
see from Eq. (|18|) that the frequency-dependent current noise is directly tied 
to the time-dependence of the second moment of the tunneled charge distri- 
bution. At long times, ((m(t) 2 )) — > flt/e, where / is the Fano factor; this 
leads to the usual expression S aa (iv = 0) = 2efl. A frequency-dependent 
S(uj) indicates that ((m(t) 2 )) deviates from this linear in time behaviour at 
short times. We thus see that in addition to being useful for calculations, 
Eq. I|18|) provides a straightforward way to interpret frequency-dependent 
current noise in terms of counting statistics. Note that unlike approaches 
using the quantum regression theorem |14l |I] , one does not need to add an 
ad-hoc term to Eq. (|18|) to obtain the correct noise in the high-frequency 
limit. 
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Figure 2. Frequency dependent current noise for the DJQP process at Ej = 10F; the 
large peak at ui = Ej is the result of oscillations in the asymmetry of the tunneled 
charge distribution. Inset: Distribution of tunneled charge p(m) for the same parameters 
at Tt = 67; the asymmetry which favours multiples of 3 is clearly visible. 



We apply Eq. (|T%|) to the DJQP process in the limit Ej » T, where 
the counting statistics are well-defined and phase-independent. As noted in 
the previous section, the coherent Josephson oscillations at frequency Ej/H 
modulates the asymmetry in p(m, t) which favours multiples of 3. This in 
turn modulates the second moment of the distribution (i.e. the greater the 
asymmetry, the smaller the width); one finds from Eq. (|14|) : 



«m 2 (t))> 



9Tt 9 
+ — 



-Ft 



cos(Ejt) / _ r4/4 



-3rt/4 



(19) 



Eq. (fT%)) then implies that there should be a peak in Sll = Srr at u = 
Ej of order e 2 Ej/T, which is much larger than the zero frequency noise 
(9/4)e 2 r. Shown in Fig. 2 are results obtained for Si(u) for Ej = 10r, 
Cl = Cr; a sharp peak is indeed visible. Similar results for the single JQP 
process were found in Ref. [I]. We see here that the peak in the current 
noise is directly related to the modulation in time of the asymmetry of the 
counting statistics by coherent Josephson oscillations. 
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5. The Intermediate Regime Ej ~ T 

The noise properties in the regime Ej ~ V are especially interesting, as 
it is in this regime where the DJQP process can be used to make a near 
quantum-limited measurement 0. In this regime, the off-diagonal elements 
of the stationary density matrix po are by no means small, and thus the 
coherent, phase-dependent contribution to the counting statistics will play 
a role. While this in itself is not a problem, the resulting status of Si(ui) 
as calculated using Eq. (fT%)) becomes unclear, as the quantity dt{(m 2 (t))) 
is phase dependent. The conclusion would thus appear to be that for res- 
onant Cooper-pair tunneling, the finite-frequency current noise is phase 
dependent, and consequently is sensitive to details of the measurement. 

We close by pointing out that the value of the zero frequency noise 
remains phase independent and unambiguous regardless of the ratio Ej/T. 
At zero frequency, each of S a p are equal; thus, we only need to know 
dt((m 2 (t))) in the large time limit. As discussed in Ref. ^To\, this is de- 
termined completely by the lifting of the zero eigenvalue of the evolution 
matrix M(k = 0) in Eq. (j3J) by non-zero k; the contribution to ((m 2 (t))) 
from other eigenmodes of M are damped away in the long time limit. 
Restricting ourselves to only this "lifted" zero-mode, we have: 

p(k,t)^e x °^ t (l + A(k,5)), (20) 

where A(k, 5) contains all fc-dependent terms involving the initial density 
matrix, and the ^-dependent eigenvectors (left and right) corresponding to 
Ao(fc); A vanishes for k — * 0. As: 

((m\t))) = ~]ogp(k,t) (21) 



it follows from Eq. (fT%|) that 

22 



d 2 

S/(0) = S aa = 2e 2 ( Ao(fc) 



k=0 



The phase-dependent term A is thus explicitly seen to play no role. Ana- 
lyzing the eigenvalue Xo(k) of M(k) in Eq. (jHJ), one find for the Fano factor 

/Eh 



J 2 



6E 2 T 2 

1 



{T 2 + 2E 2 Y 



(23) 



In both the limits Ej <C Tor Ej ^> T , / — » 3/2. This can be understood as 
an effective charge- each Cooper-pair plus quasiparticle transition transfers 
on average 3/2 electrons per junction. Interestingly, when Ej ~ T, f drops 
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below 1/2, reaching a minimum of 3/8 when Ej = T/\/2. This behaviour 
is reminiscent of double tunnel junction systems, or of a normal SET; in 
both cases, there are two rates involved in transport, and the Fano factor 
reaches a minimum of 1/2 when these rates are equal. The behaviour in such 
systems can be understood completely classically by the self-averaging that 
occurs when one sequentially combines two independent Poisson processes. 
In contrast, the fact that / drops below 1/2 for the DJQP process is a direct 
consequence of the coherence of Cooper-pair tunneling. To underscore this 
point, one can calculate / using the incoherent model of Eq. ©. While 
this procedure yields the correct value of (J) and the correct value of / for 
extreme values of Ej/T, it only gives a minimum / of 1/2. One finds: 



/ incoherent 2 



2#2r 2 

1 J — 



(l* + 2E°) a 



f+ ^ (24) 



Note that the incoherent calculation of / is always an overestimate; as 
expected, coherence between charge states suppresses the zero-frequency 
current noise, as it tends to make tunnelling events more regular. 



6. Conclusions 



We have studied the counting statistics and finite-frequency shot noise 
of the DJQP process in a superconducting single electron transistor, at- 
tempting to clarify some of the remarkable noise features found in previous 
studies. In general, the counting statistics are phase-dependent, despite 
the fact that the average current has no phase dependent. An interesting 
question remains how this phase dependence impacts the finite frequency 
noise, given Eq. (fTH|) which directly relates the two quantities. 
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